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Abstract 

The pumping intensity (/) dependence of the photoluminescence (PL) spectra of perfectly lat- 
erally two-dimensionally ordered SiGe quantum dots on Si(OOl) substrates was studied. The PL 
results from recombinations of holes localized in the SiGe quantum dots and electrons localized 
due to the strain field in the surrounding Si matrix. The analysis of the spectra revealed several 
distinct bands, attributed to phonon-assisted recombination and no-phonon recombination of the 
excitonic ground state and of the excited excitonic states, which all exhibit a linear I dependence 
of the PL intensity. At approximately 7>3Wcm -2 , additional bands with a nearly quadratic / 
dependence appear in the PL spectra, resulting from biexcitonic transitions. These emerging PL 
contributions shift the composite no-phonon PL band of the SiGe quantum dots to higher energies. 
The experimentally obtained energies of the no-phonon transitions are in good agreement with 
the exciton and biexciton energies calculated using the envelope function approximation and the 
configuration interaction method. 

PACS numbers: 73.21.La,73.22.Lp,78.55.Ap,78.67.Hc 
Keywords: quantum dots, SiGe, biexciton, k ■ p 
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I. INTRODUCTION 



SiGe quantum dots (QDs) are of interest, since they provide a promising way towards an 
infra-red light source operating on telecom wavelengths, integrable into the present Si-based 
technology- Although these structures have been studied for quite a long time,- 1 ^ there are 
only a few reports on the excitonic structure of their photoluminescence (PL) spectra." The 
large size of the QDs and an inhomogeneous Ge distribution due to intermixing and clustering 
lead to a pronounced broadening of the PL lines of QD ensembles.- Additional broadening 
arises from the indirect nature of the optical transitions both in real and reciprocal space. 
Due to the broad spectra, an unambiguous assignment of the electronic QD transitions of 
randomly nucleated QDs grown on planar Si (001) has not been reported so far. For the 
same reasons, up to now only continuous linear and sub-linear shifts of the PL emission 
energies with increasing excitation intensity (J) could be observed both for SiGe QDs&2 and 
quantum wells (QWs),— ~— which were attributed to state filling effects and carrier induced 
band-bending as a consequence of the spatial separation of the excited electrons and holes. 

It has been shown that by a precise positioning of the QDs on pre-structured Si sub- 
strates,— a significant narrowing of the PL spectra can be achieved,-^— making a more 
detailed analysis of the dependence of the PL spectra on the excitation intensity possible. 

However, no decomposition of the PL spectra in terms of individual excitonic recom- 
binations has been reported so far for laterally ordered SiGe QDs. Here we provide an 
identification of both bound excitonic and biexcitonic transitions in the PL spectra of en- 
sembles of ordered SiGe QDs, using the linear and quadratic excitation intensity dependence 
of those transitions. This identification is facilitated by the clear observation of additional 
PL emission lines that appear with increasing / at higher photon energies as opposed to the 
continuous shift of the PL spectra reported up to now. The assignments are supported by 
the results of exciton and biexciton energy level calculations. 

Extensive work on the exciton-exciton interaction energies has been reported for self- 
assembled epitaxial III-V compound quantum dots^~— as well as for chemically synthesized 
colloidal nano crystals^— for material systems with both type-I and II band alignment. The 
understanding of the multi-exciton interaction and the exciton relaxation dynamics has been 
shown to be important for the realization of hot carrier extraction and optical gain.— >^ 

In this work the exciton-exciton interaction as well as the exciton excited state energies 
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are determined experimentally for SiGe QDs epitaxially embedded in a Si matrix. In this 
material system, the holes are confined in the SiGe QDs and the electrons to the tensile 
strained regions in the Si matrix surrounding them, forming a type-II system. Our results 
show that for the ordered, highly uniform SiGe QDs investigated in this work, the energy 
splitting between the biexciton and exciton emission line is larger than the line width, which 
is a prerequisite to achieve optical gain for single exciton states.— 

This paper is organized as follows: in SecHU experimental details on SiGe island growth 
and on the PL measurements are presented. In Sec. lIIIt the analysis of the PL spectra, 
their excitation intensity dependence, electronic structure calculations and the resulting 
assignment of the exciton, excited exciton, biexciton and excited biexciton energies are 
given. In Sec.lIVt the conclusions are presented. 

II. EXPERIMENTAL 

Two samples with different average Ge content in the SiGe QDs were investigated (SI 
and S2). High resistivity p-type Si (001) substrates were pit patterned by nanoimprint 
lithography and subsequent pattern transfer into the substrate by reactive ion etching. The 
pattern periods were 300 nm for SI and 170 nm for S2, and the patterned areas were 3x3mm 2 
for both samples. The resulting pits had diameters and depths of 160 nm and 47 nm for 
SI, and 120 nm and 35 nm for S2, respectively. After standard cleaning, the wafers were 
in-situ degassed in the solid source molecular beam epitaxy chamber for 40 min at 720 °C. 
Hereafter, a 45 nm thick Si buffer layer (growth rate R = 0.6 A/s) was deposited at a 
substrate temperature that was increased from 450 °C to 550 °C followed by the deposition 
of 6 monolayers (ML) of Ge at 690 °C (R = 0.05 A/s) for SI and 8.3 ML of Ge at 625 °C 
(R = 0.025 A/s) for S2. Subsequently, a 50 nm or a 10 nm thick Si capping layer was 
deposited at a temperature as low as 300 °C in order to avoid unwanted Si incorporation 
and QD shape transformations for SI or S2, respectively— Atomic force micrographs of the 
ordered and capped islands are shown in Fig. [T]for SI and S2. 

In the PL measurements, performed at 4.2 K, the samples were excited by an Ar + laser 
tuned to the wavelength of 457.9 nm, focused by a lens to a circular area with 400 fim 
diameter (~ 10 6 irradiated QDs). The excitation intensities ranged from 0.15 Wcm" 2 to 
7.90 Wcm" 2 for sample SI and from 0.15 Wcm" 2 to 3.95 Wcm -2 for sample S2. After being 
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dispersed by a grating monochromator, the PL light was detected using an InGaAs line 
detector at the temperature of —100 °C. 

III. RESULTS AND DISCUSSION 

A. Analysis of the photoluminescence spectra 

In Fig. E£a)— (c) the PL spectra of SI are shown for three excitation intensities (J = 0.25, 
0.49, 4.94 Wcm -2 ) together with their decomposition into various lines according to the 
fitting procedure described below. It is evident that with increasing / additional lines 
appear in the high-energy shoulder of the PL spectra. The observed behavior can not 
be described by a continuous shift of the emission energy caused by band-bending as a 
consequence of electron-hole separation due to the type-II band alignment in the SiGe QDs. 
For the complete range of excitation intensities the PL spectra of SI and S2 are shown 
on a logarithmic scale in Figure |3J Due to the lower growth temperature for the islands 
in sample S2 as compared to SI, both the average and the maximum Ge content in S2 is 
larger.- Thus, the island related PL spectrum of S2 is observed at lower energies than for SI. 
For both samples, the appearance of additional lines with increasing / (as shown in detail in 
Fig. [2]) is evident over an excitation intensity range between 0.15 and 7.9 (3.95) Wcm -2 for 
SI (S2). In the following, a quantitative description of our observations in terms of excitonic 
contributions will be given. 

In order to identify the transitions contributing to the individual PL spectra as a function 
of /, the spectra were fitted using the Gauss-Lorentz (GL) profiles, employing the method 
of the rational approximants.— Three parameters were used to describe each contribution: 
the resonant energies (E ^), the oscillator strengths (Fj) and the Gaussian widths (T Gti ), 
where i indexes the GL profiles. For the spectral region of phonon replicas below 885 meV 
(828 meV) for SI (S2) (see also Fig. [2]), the values of all these parameters were adjusted to 
fit the data, while for the region of no-phonon transitions the values of E Q ^ were fixed and 
only those of Fi and Tc,i allowed to vary in the fitting routine. The fixing of the values of 
i?o,i has been motivated by our experimental finding that in the no-phonon region the lines 
indeed appear at E a i and are fitted more precisely with increasing I. This is in agreement 
with the assumption that every profile in the region of no-phonon transitions corresponds 
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to one excitonic or biexcitonic transition. The width of the Lorentzian contributions T i,i to 
the fitted line width was fixed at the small value of 0.001 meV for all the fitted GL profiles, 
i.e., the profiles were treated as almost purely Gaussian ones. The Lorentzian width is 
negligible at low temperatures because the dominant spectral broadening is inhomogeneous, 
originating mainly in the statistical variation of the QD structure. 

We have fitted every PL spectrum using the smallest number of profiles needed for a 
reasonably good fit, assessed by the residual sum of squares and the correlation coefficients 
between the fitted parameters. The values of the parameters obtained for I n were used 
as starting values for fitting the spectrum measured at I n +i. Here I n and J n+ i denote 
subsequent excitation intensities in the series shown in Fig. [3j If, for the given number of 
profiles the best agreement at I n+1 was significantly worse than that at I n , an additional 
profile was added. In this manner we have identified, for both samples, the number of profiles 
of the phonon assisted part of the spectra to be 3 or 4, and the number of the no-phonon 
part ranging from 1 to 3. As an example, we show in Fig. [2] three selected PL spectra of the 
sample SI along with the decomposition into the individual GL profiles. The values of E 0ti 
are summarized in Tab. HI 

The phonon assisted transitions were attributed to the various SiGe phonon modes typi- 
cally observed in the PL spectra of bulk SiGe samples^ and quantum wells.— With increasing 
value of Eq these are the Si-Ge transverse optical (TO), the Ge-Ge TO (for SI resolved only 
at / > 1.48 Wcm" 1 ), the longitudinal acoustic (LA) and the transverse acoustic (TA) (for 
S2 resolved only at / > 0.20 Wcm" 1 ) phonon assisted transitions. The Si-Si TO phonon 
assisted transition was not identified in our spectra; its contribution may overlap with the 
band attributed to the Si-Ge TO phonon replica. With respect to the no-phonon transitions, 
for the sample SI (S2) and for / in the range from 0.15 to 0.40 Wcm -2 (0.15 to 0.30 Wcm -2 ) 
a single profile with a fixed value of E = 888 meV (Eq = 832 meV) was used. For / in the 
range from 0.49 to 0.62 Wcm -2 (0.35 to 0.99 Wcm -2 ), a second profile with E = 893 meV 
(Eq = 842 meV) was added and another profile with Eq = 898 meV (Eq = 848 meV) was 
added for / in the range from 0.99 to 7.90 Wcm -2 (1.48 to 3.95 Wcm -2 ). The purely 
electronic (no-phonon) transitions are attributed to the ground state exciton (Xo), excited 
exciton (Xi), and ground state biexciton (XXo) in the case of SI, see Fig. [2} and to Xo, 
XX and excited biexciton state (XXx) in the case of S2. This interpretation is based on the 
observed / dependence of the oscillator strengths Fj (see subsection IIIIBI) . and the results 
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of electronic structure calculations (see subsection IIII Cj) . 



B. Excitation intensity dependence 

The I dependencies of the oscillator strengths and of the Gaussian full widths at half 
maximum (FWHM ) of the fitted profiles of the no-phonon transitions are displayed in Fig. [5] 
and its insets, respectively. In order to suppress the uncertainties of the pumping intensities 
all Fi values shown in Fig. [4] for sample SI (shown in panel a) [S2 (panel b)] are normalized 
to the oscillator strength of the 888 meV (832 meV) band, which corresponds to the lowest 
excitonic state (Xq, with Fx linear in J). Figure Ufa) reveals an approximately linear 
(~ ji-03-j dependence of F of the 893 meV band and a superlinear (~ J 2 - 09 ) dependence of 
the 898 meV band of the sample SI. For the sample S2, both the 842 meV band and the 
848 meV one exhibit a superlinear, (~ J 170 ) and (~ J 192 ) dependence of F, respectively, 
as shown in Fig.|4^b) (Note, that due to the linearity of the normalization factor Fx with 
respect to I, quadratic (linear) powers of / appear as linear (constant) functions in Fig. [4]). 

Figure H] also shows the I dependence of the sum of the oscillator strengths of the phonon 
replicas. This dependence is also slightly superlinear, ~ 7 111 and ~ I 1A2 for SI and S2, 
respectively. We interpret this finding as an indication for contributions of phonon replica 
lines of non-linear exciton emissions. Note, that for the sample S2 the relative magnitude of 
the contribution to the phonon-replica bands superlinear in / is larger than for the sample 
SI. This may be due to the fact that the ratio of the sum of F for the no-phonon bands 
having a superlinear I dependence to that of those having a linear I dependence is higher for 
S2 than for SI. The broad features in the I dependence of F and FWHM (the latter shown 
in the inset of Fig. H]) are fitting artifacts caused by high correlations of the parameters 
Fi and T^ j, respectively. Note, that the values of the FWHM of the no-phonon bands are 
comparable to the FWHM of the X band observed for SiGe bulk crystals (8 meV, see Fig. 10 



of Ref. l29l ). However, they are considerably larger for S2 than for SI (by a factor of ~ 2). 

The insets of Fig. H] show that the differences between the various exciton resonance ener- 
gies are comparable to the widths of the GL profiles, i.e. despite the excellent homogeneity 
of the QD ensemble these differences are just beyond the experimental resolution limit. 
Thus, any additional inhomogeneous line broadening, such as that occurring in randomly 
nucleated islands, hindered the discrimination of various excitonic contributions to the PL 



spectra in previous studies. 



C. Electronic structure calculations 

To get a better insight into the origin of the no-phonon transitions, a series of calculations 
of the electronic structure was performed using the following two-step approach. First, the 
single particle wavefunctions were obtained by the nextnano++ solver— using the single 
band effective Schrodinger equation for the electron states in the A valley of the lowest 
conduction band and the six-band envelope function method for the hole states. The electron 
and hole states were thus treated as decoupled. This approximation is justified considering 
the energy difference of ~ 750 meV between the extrema of the confinement potentials of 
electrons and holes. Second, these calculated wavefunctions were used to construct a basis 
set for the configuration interaction (CI) calculations.— The evaluation of the Coulomb 
matrix elements is facilitated by the orthogonality of the periodic parts of the Bloch waves 
of the bottom of the conduction band in Si and the top of the valence band in both Si and 
Ge.^ 

Next we describe the model structures. The SiGe QDs of sample SI were defined on the 
rectangular grid and approximated by cones with base diameter of 122.8 nm and a height of 
15.3 nm as sketched in Fig. [5j For sample S2, a similar structure was used with slightly dif- 
ferent dimensions of 118 nm and 15 nm deduced from the AFM measurements on uncapped 
islands. In the model structure representing the sample SI (S2) the Ge content linearly 
increases from 0.277 (0.34) at the QD base to 0.43 (0.49) at its apex; this is motivated by 
results of Ref. Q. The Ge concentration profiles were chosen to warrant agreement between 
the measured and calculated values of the transition energy of the lowest no-phonon transi- 
tion. For more information on the model structures, see Fig. [51 for the material parameters, 
see Tabs. [TTJ IIHI We have found that the uncertainty in the values of the deformation poten- 



tials in Si and Ge (estimated from the spread of the values published in Refs. I34H37I) leads 
to an uncertainty of ~ 20meV in the calculated value of the energy of X . However, the 
differences between the energy of X and those of the other excitonic complexes are almost 
independent of the choice of the deformation potential parameters. The spacing of the grid 
used in the calculations was set to 4 nm in both lateral and vertical directions except for 
a cuboid around the QD apex (shown in Fig. [5] by the broken line) where the spacing was 
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0.5 nm in all directions. The Schrodinger equation was solved only in this subspace owing to 
the expected positions of the electron and hole states.— On the other hand, the minimization 
of the strain energy was performed in the whole simulation space. For both calculations von 
Neumann boundary conditions were used in both lateral and vertical directions. 

D. Assignment of the PL bands 

The calculated values of the transition energies of the states Xo, Xi and XXo for the 
sample SI are indicated in Fig. [2] by the broken and full vertical lines along with the prob- 
ability densities of the single electron wavefunctions [inset of panel c)]. Note that for both 
samples the hole wavefunctions were predominantly composed of heavy holes (96%); the 
electron wavefunctions, from which the lowest excitonic (and biexcitonic) states were com- 
posed, belonged purely to the lower lying A z conduction band valleys oriented with their 
main axis along the [001] growth direction. This is due to the difference of about 18meV 
between the energies of the lowest A z and the four A xy single electron states and almost no 
spatial overlap between the A xy state and the hole states (these results have been obtained 
by calculations involving the whole simulation space with a less dense grid).— 

The resulting energy of the excitonic ground state (X ) for the sample SI, calculated 
using the CI method with 6 electron and 4 hole basis states, was found to be 887 meV 
(the corresponding single particle energy difference between the electron and hole states is 
896 meV), close to the observed value of 888 meV, and that of the first exited excitonic state 
(Xi) was found to be 894 meV. The biexciton (XXo) was found to be shifted to higher energies 
with respect to the exciton ground state by 11 meV, a value which is in very good agreement 
with the experimentally observed energy difference of 10 meV between the first and the third 
no-phonon GL profile having the linear and a superlinear / dependence of F, respectively. 
Also the observed difference of 5 meV between the energies of the first and the second no- 
phonon profile is in reasonable agreement with the calculated value of £'(X 1 ) — E(K ) of 
7meV. Thus, we assign the 888 meV and 893 meV profiles to the recombination of Xo and 
Xi, respectively, and the 898 meV profile to the recombination of XXo- This assignment is 
corroborated by the observed dependence of the strength of the respective GL profiles on I 
as discussed in section [Til Bl 

For the sample S2, the energy of X , calculated using the CI method with 6 electron 
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and 4 hole basis states, was found to be 837 meV (the corresponding single particle energy 
difference between the electron and hole states is 847 meV), close to the observed value of 
832 meV. The calculated values of the blueshift of the biexciton ground state XX and of the 
biexciton excited state XKi are 12 meV and 18 meV, respectively, in good agreement with the 
experimentally observed energy differences of 10 meV and 16 meV between the second and 
the first no-phonon band, and between the third and the first no-phonon band, respectively. 
Thus, we attribute the 832 meV, 842 meV and 848 meV bands to the recombination of 
Xo, XXo and XXi, respectively, again in agreement with the observed dependence of their 
strength on / as discussed in section [III Bl The excited excitonic state Xi observed in sample 
SI was not identified here. This is probably due to the larger inhomogeneous broadening of 
the GL profiles in S2 compared to SI. The Xi band might contribute to the second no-phonon 
profile, causing a slightly lower magnitude of its superlinear / dependence seen in Fig. H^b) . 
Also, excited state surface trapping, as observed in Ref. |24j for colloidal nanocrystals, might 
quench more efficiently the Xi emission in S2 as compared to SI due to the thinner capping 
layer of S2. The assignment of the 842 meV and 848 meV bands to XXo and XXi is not 
unique as a similarly satisfying agreement with the experimental data may be achieved by 
attributing these to the recombination of the groundstate trion and a XX at an even higher 
energy, respectively. However, in our calculations this would require the assumption of a Ge 
concentration profile in the QD, with an - for the ordered SiGe QDs - unrealistically large 
Ge accumulation of up to 60-70% at its apex. Evidence for such a large Ge concentration 
in the apex of SiGe transition dome and dome islands have been found up to now only in 
randomly nucleated islands grown on planar substrates, but not in ordered ones.- 1 ^ 1 ^ 

Our calculations and experiments show that the exciton-exciton interaction is pro- 
nouncedly anti-binding in SiGe QDs, resulting in ~ 10 meV higher biexciton transition 
energies as compared to the exciton ones. A similar antibinding exciton-exciton interac- 
tion was invoked in the interpretation of Ge hut-cluster absorption spectra.- Since also in 
InAs/GaAsSb^, GaSb/GaAs^ and in InP/GaAs2i quantum dot systems with type-II band 
alignment a large anti-binding exciton-exciton interaction was found, we conclude that the 
antibinding character of the exciton-exciton interaction is characteristic of type-II systems 
with spatially separated electrons and holes.— Our conclusion is also supported by ex- 
tensive work on exciton-exciton repulsion in colloidal semiconductor nanocrystals reported 



in Ref. 



23 



and 



26 
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IV. CONCLUSIONS 



In conclusion, we have performed extensive intensity dependent PL measurements on 
two ensembles of two dimensionally ordered SiGe QDs, whose excellent homogeneity is 
due to the controlled growth on prepatterned Si substrates. The spectra were decomposed 
into a series of distinct bands with characteristic excitation intensity dependencies of their 
oscillator strengths. Electronic structure calculations were performed using the nextnano++ 
solver and the calculated wavefunctions were used as a basis set for configuration interaction 
calculations. Based on these calculations, the transition energies of the Xo,Xi,XXq and 
XX i states were compared with the experimentally observed PL bands. Together with 
the experimentally observed excitation intensity dependence of the various PL bands, the 
excitonic and biexcitonic recombinations are identified in this type-II quantum dot system. 
To the best of our knowledge this is the first clear evidence for the formation of excitonic 
complexes in this system, favored by the zero-dimensional nature of the SiGe/Si QDs. The 
determined exciton interaction and excited state energies are important for the application 
of SiGe QDs in Si photonics. 
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FIG. 1. Atomic force micrographs of the samples (a) SI (pit period 300 nm) and (b) S2 (pit 
period 170 nm). The scanned area was 1.5 x 1.5 /mi for both samples. The average height of the 
QDs in sample SI (S2) is 15.3 (15) nm and the QD base diameter is 122.8 (118) nm. 
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FIG. 2. Measured (open squares) and fitted (red curve) PL spectrum of the sample SI for (a) 
/ = 0.25 Wcnr 2 , (b) I = 0.49 Wcm~ 2 and (c) I = 4.94 Wcrn" 2 . The residual error (blue curve) 
has been multiplied by a factor of 2. The individual GL profiles are attributed to phonon replicas 
(dotted curves), excitonic (Xq, X\, broken curves) and biexcitonic (XXq, full curve) states. The 
calculated values of the energies of the excitonic (biexcitonic) states for the model dome-shaped 
dot (see Fig. ED blueshifted by 1 meV are displayed by dashed (solid) vertical lines. The inset 
of panel (c) shows the calculated probability densities (contours of \^ 2 \ = 0.1 nm -3 ) of the hole 
ground state (blue), A xy (green) and A z (red) electron ground states and their location within the 
dot (its surface is represented by the light green surface) as obtained by the nextnano++ simulation 
suited 
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FIG. 3. Measured PL spectra (the values of the pumping intensity are given next to the right 
vertical axis) for samples (a) SI and (b) S2. Dashed (solid) vertical lines correspond to the resonance 
energies Eq of the fitted GL profiles attributed to excitons (biexcitons). The dotted vertical line 
in b) corresponds to the detector cutoff energy of 780 meV. 
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FIG. 4. Excitation intensity (/) dependence of the oscillator strength (Fi) normalized to the 
respective squared resonance energies Eq^ and divided by the oscillator strength Fx of the no- 
phonon transition with the lowest value of Eq, i.e. of the Xq profile, for the no-phonon transition 
GL profiles in the sample SI (a) and in the sample S2 (b). The sums of F^ of the phonon replicas 
are displayed by open triangles. The thick lines represent fits of the measured data to the linear 
functions a\ ■ I 0,2 and the thin ones are guides to the eye. Note, that due to the linearity of 
the normalization factor Fx with respect to /, quadratic (linear) powers of / appear as linear 
(constant) functions. The graphs in the insets show the / dependence of the FWHM of the fitted 
Gaussian lines. The same symbols as in the main diagram are used except for the dependencies of 
the profiles with the lowest value of Eq which are displayed by open circles. 
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FIG. 5. A cut through the simulation space for the SiGe QD of sample SI (S2), containing its 
vertical axis. The QDs were approximated by cylindrically symmetric cones with the dimensions 
shown in the figure. The Ge content inside the QD varies linearly from 0.277 (0.34) at the base of 
the cone to 0.43 (0.49) at its apex. The dotted rectangle denotes the space where the Schrodinger 
equation was solved and the grid spacing was set to 0.5 nm. 



19 



TABLE I. Energies of the fitted GL profiles and calculated energies of the no-phonon excitons for 
samples SI and S2 along with the energies of the phonon replicas taken from Fig. 8 of Ref. [29. 



The energies of the phonon replicas and the excitonic complexes are given relative to the corre- 
sponding value of Xo. The displayed values of the energies of the phonon replicas were taken from 
the corresponding GL profiles of the spectra for the lowest excitation intensities. The estimated 
uncertainties of the energies of the fitted GL bands are ~ 3 meV for the no-phonon bands and 
~ 6 meV for the phonon assisted ones. For the uncertainties of the calculated values of Xq see 
text, the uncertainties of the calculated energies of the other excitonic complexes is ~ 2 meV. The 
uncertainty of the phonon energies derived from Fig. 8 of Ref. Q is ~ 1 meV. 





X 


Xi XXq XXi 


TO Si-Ge TO Ge-Ge LA TA 


SI fit (meV) 
SI theory (meV) 


888 ±3 
887 ± 20 


+5 ±3 +10 ±3 
+7 ±2 +11 ±2 


-52 ±6 -39 + 6 -26 + 6 -9 + 6 


S2 fit (meV) 
S2 theory (meV) 


832 ±3 
837 ± 20 


+10 + 3 +16 + 3 
+8 + 2 +12 + 2 +18 + 2 


-46 ± 6 -34 ± 6 -19 ±6-9 + 6 


Ref. 29 (meV) 






-49 -35 -30 -10 
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TABLE II. Description of the material parameters used in the calculations, 71, 72 and 73 are the 

Luttinger parameters.— 

Parameter Description 

a lattice constant 

a exp lattice thermal expansion coefficient 

Cn elastic constant 

C12 elastic constant 

C44 elastic constant 

e r static dielectric constant 

A-valley longitudinal electron effective mass 
TOAt A-valley transversal electron effective mass 

Eq bandgap 
a Varshni parameter 

(3 Varshni parameter 

E v valence band offset 

Ao spin-orbit split-off energy 

a c absolute deformation potential for conduction band A-valley 

a u uniaxial shear deformation potential of the conduction band A-valley 

a v absolute deformation potential for valence band 

a n fe uniaxial shear deformation potential b of the valence bands 

a u d uniaxial shear deformation potential d of the valence bands 

L Dresselhaus parameter;— L = —71 — 472 — 1 

M Dresselhaus parameter^ M = 272 — 71 — I 

iV Dresselhaus parameter^ TV" = —673 
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TABLE III. Values of the material parameters used in the calculations. The unit mn represents 
the free electron mass. 



Parameter Unit 


Si 


Ge 


Sii_ x Ge x 


a 


A 


5.4304 [43] 


5.6579 [43J 


linear 




A/K 


1.8138xl0~ 5 [43] 


5.8xl0~ 5 [44] 


linear 


Cn 


GPa 


165.77 [43] 


128.53 [43] 


linear 


C12 


GPa 


63.93 [43J 


48.26 [43] 


linear 


( j A A 


GPa 


79.62 [43J 


66.80 [43] 


linear 






11.7 [45] 


16.2 [43] 


linear 


7 77 a 7 


777 n 


0.916 [45] 


1-350 [21] 


linear 


777, A + 


771 n 


0.190 [45] 


0.290 [21] 


linear 


En 


eV 


1.17 [46] 


0.931 [29] 


0.931x + 1.17(l-x) - 0.206x(l-x) 




eV/K 


0.473xl0" 3 [44] 


0.4774xl0~ 3 [31] 


linear 


3 


K 


636 [44] 


235 [2U 


linear 




eV 


1.090 [37] 


1.67 [37] 


linear 


An 


eV 


0.044 [43] 


0.30 [43J 


linear 


1 1 ■(... 


eV 


3.40 [37] 


0.14 [37] 


linear 




ev 


9.16 [35] 


9.42 [35] 


linear 


a v 


eV 


2.05 [37] 


-0.35 [37] 


linear 


a u b 


eV 


-2.10 [47] 


-2.86 [48] 


linear 


dud 


eV 


-4.85 [47] 


-5.28 [48] 


linear 


L 


h 2 /2m 


-6.69 [49J 


-31.34 [50] 


linear 


M 


h 2 /2m 


-4.62 [49J 


-5.90 [50] 


linear 


N 


h 2 /2m 


-8.56 [49J 


-34.14 [50] 


linear 



[29] 
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